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The efficiency, robustness and reliability of recent numerical methods for finding solutions to flow problems have 


given rise to the implementation of computational fluid dynamics (CFD) as a broadly used analysis method for 


engineering problems like membrane separation system. The CFD modeling in this study observes steady and un- 


steady (transient) heat flux and temperature profiles in a polymeric (cellulose acetate) membrane. This study is 


novel due to the implementation of user defined scalar (UDS) diffusion equation by using user-defined functions 
(UDFs) infinite volume method (FVM). Some details of the FVM used by the solver are carefully discussed when 


implementing terms in the governing equation and boundary conditions (BC). The contours of temperature due to 


high-temperature gradient are reported for steady and unsteady problems. 
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Introduction 


Computational methods have the ability for improving 
the understanding of flow problems in membrane separa- 
tion systems. They have the proficiency of giving data on 
flow conditions at any point of the geometry without 
troubling the flow. Furthermore, the implementation of 
mathematical modeling expressively decreases the costs, 
risks, and time-related to performing repetitive experi- 
ments. One numerical method implemented for simulat- 
ing fluid flow is known as computational fluid dynamics 
(CFD) [1-4]. CFD has become a more commonly im- 
plemented method in the research area of membrane 
science [5, 6] Several researchers are using this tool in 
order to obtain understanding into the phenomena taking 
place inside membrane modules, to support the design 
process and increase the performance of modules. As an 
analysis method, CFD gives the permission to change 
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fluid properties, operating conditions and geometric cha- 
racteristics of the flow channels in a flexible but defined 
way. The geometric parameters of the channels simulated 
can be changed deprived of the requirement to make and 
check new meshes. This signifies a major benefit of the 
CFD tool over traditional experimental methods. Also, 
flow data can be described at any point and time through- 
out the simulation, and the assessment of these flow va- 
riables is conceivable deprived of any trouble to the 
original flow [7, 8]. Because of its easiness and minor 
computational requirements, many scholars have chosen 
CFD for modeling membrane in two and three dimen- 
sions, as is showed by the quantity of effort carried out in 
this area [9-13]. Consequently, CFD is currently reco- 
gnized as a dependable technique, and progressively un- 
steady flow studies are developing, mostly finding decent 
similarity within calculated and experimental results [14]. 

Organic and inorganic materials are used to synthesize 


membranes for gas separation. In late 1970s commercial 
polymeric membranes i.e. polydimethylsiloxane, cellu- 
lose acetate, polyether sulphone etc. for gas applications 
became available [15]. Membranes based on polymers 
establish extensive use in the separation of N2, O2, and 
CO, at lower temperatures (<100°C). Though, several 
engineering methods include the production of gases at 
high temperatures. One case is the production of syngas 
(a mixture of carbon dioxide and hydrogen) through coal 
gasification or natural gas reforming [16]. Preliminary 
research directed to different ideas on the consequence of 
gas polarization which is linked with the reduction of the 
additional permeable gas species through the membrane. 
This consequence allows the increase of the less permea- 
ble gas on the feed side and have a tendency to increase 
in significance for membranes providing very high selec- 
tivities [17]. But, it is reported that at high temperatures 
the gas-through-gas diffusivity in the feed stream is quite 
a few orders of higher magnitude than the diffusion of 
gases through the membrane, therefore decreasing con- 
centration polarization [18]. 

Membrane gas system is an important separation 
process of chemical engineering. In this research, we 
focused on heat transfer through the domain of a 2D rec- 
tangular membrane using general scalar diffusion equa- 
tion with mixed boundary conditions. This research cov- 
ers the solution of steady and unsteady state problems 
with the formulation of boundary conditions for first or- 
der discretization technique. The UDF is integrated with 
FVM based CFD solver (ANSYS FLUENT™) for the 
solution. The contours of rectangular membrane domain 
show the behavior of a temperature through a polymeric 
(cellulose acetate) membrane. The heat transfer prob- 
lem of flow in a membrane is a less focused area, espe- 


cially for high-temperature gradient due to its complexity. 


This study is an effort to solve user defined scalar (UDS) 
diffusion equation to solve complex thermal diffusion 
phenomena in a membrane. 


CFD modeling 


The heat transfer modeling and the phenomena of 
transport mechanisms are permanently linked with the 
fundamental properties of the material used to synthesize 
a membrane. A common principle of membrane opera- 
tion is that there must be a driving force for a gas mole- 
cule to diffuse from the feed domain to the permeate do- 
main. This driving force is generally expressed as a tem- 
perature gradient in the case of heat transfer in polymer- 
based membranes. The membrane is described as a ze- 
ro-thickness wall. The first cell in the first row of a mesh 
on the feed side is considered as a source. Introducing the 
negative source and source term at two sides of the 
membrane permits the gas disappear from the feed do- 


main via the membrane surface and emerge at the per- 
meate domain. The mesh of membrane is developed by 
using GAMBIT ™. Grid size analysis is carried out using 
different mesh intervals. All the simulation results did not 
show any major difference. The selected membrane mesh 
is shown in Figure 1. The mesh consists of 800 cells with 
the dimension of 0.05 m x0.1 m. 
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Fig.1 2D mesh of a membrane 


FVM solver (ANSYS FLUENT™) solves the trans- 
port equation for a user-defined scalar (UDS) in the same 
way as it solves the transport equation for a scalar in the 
core equations, such as a species mass fraction. The UDS 
capability can be used to implement a wide range of 
physical models. In this study, a model is developed to 
solve a general scalar diffusion equation (1) with the 
possible types of boundary condition (BC) at the boun- 
dary (or a part of the boundary) of the domain [19]. 


cP_v(TV9)=5,,in2,1>0 (1) 


The possible types of boundary conditions (BC) are as 
follows: 


Dirichlet BC: Do, Neumann BC: -T (0¢/0n) = q and 
mixed BC: -I'(0¢/én)= h,(¢-¢,) 

Here, Do, qo and ø, are constant values. The boun- 
dary conditions for left-wall and top-wall are q= 
-I'(0¢/6n)=0 and ¢ = 30 respectively. The following 


equation gives mixed boundary condition for right-wall 
and bottom-wall. 


-(r2 =h,(¢-¢,) (2) 


Steady State modeling 
Consider a steady-state scalar equation with constant 
T and zero source term (Sy= 0) as follows: 
-V(I'V¢)=0 (3) 
This is the Laplace's equation. Solving the Laplace’s 
equation using the FVM UDS solver (ANSYS FLUENT™) 
does not necessarily require user-defined functions 
(UDFs). We can activate the UDS from the graphical user 
interface. FVM solver UDS provides only Dirichlet and 
Neumann conditions for the boundaries. Hence, we need 


to use the UDF to apply the mixed boundary condition 
for the UDS equation. 


Formulation of mixed boundary conditions 


For a generic cell c0 adjacent to the boundary, the dif- 
fusive flux across the boundary face f of the cell is ex- 
pressed as follows: 


-f r(2)as- fh (¢-¢, dS (4) 
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Fig.2 Schematic of a cell adjacent to a boundary 


Using the mid-point rule of surface integral, the diffu- 
sive flux can be approximated as: 


r (3$) Ar ahr to) 4s (5) 
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In FVM solver, the diffusive flux is approximated in 
two parts; in the first part primary gradient is evaluated 
implicitly along the line connecting the cell centroid c0 to 
the centroid face f. It is corrected by a secondary gradient 
(or cross diffusion) term evaluated explicitly by the gra- 
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FVM solver provides specific macros for UDF which 
define the necessary geometrical variables of the cell and 
also calculate the secondary gradient term in equation (6). 
If we designate the secondary gradient term in equation (6) 
Po and A, A/A,e, as Abe, equation (5) and (6) can be 
written as 
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¢ can be expressed as follows: 
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The mixed boundary condition for the UDS is ready to 
be specified by boundary profile grin equation (8) through 
the UDF macro [20]. 


Unsteady state modeling 


The unsteady (transient) user-defined scalar equation 
is as follows: 
0g 
ea ViIVe) =O. .t> 0 (9) 
It has to be solved with the same boundary conditions 
and given an initial condition. When the transient term 
is 0(p¢) [ôt , FVM can readily handle this unsteady UDS 


term when you enable the unsteady solver (either first or 
second-order). But if the transient term of the user's equ- 
ation is not the same as the given form, you must supply 
the unsteady term through the specific UDF macro. For 
example, in equation (9), c is a variable. For an unsteady 
heat conduction problem in dimensional form it repre- 
sents pc, of the solid material. 

In finite-volume methods, the transient term is first 
approximated by first or second order finite difference 
expression and then integrated with respect to the cell 
volume. The FVM solver expects this transient term to be 
moved to the right-hand side of the governing equation 
and included in the discretized equation as a source term. 


First-Order Unsteady Formulation 
The first-order finite-difference backwards differenc- 
ing approximation gives: 
op g'-9" 
t At 
Where, ¢’ is the value of ¢ at the current time interval, 


gn: is the value at previous time level and Aż is the time 
step size. Hence 
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where, AV is the volume of each individual cell, Apu is 
the coefficient multiplying ø” (the implicit part) and S, 
the explicit part (because it is expressed by known values 
at the previous time-step). We can rewrite the volume 
integral of the unsteady term as: 


-Í Lay = A,,6" +S, (12) 


and it is ready to be coded in UDF as the unsteady macro. 
In most existing types of CFD software based on FVM, 
this problem can be resolved by introducing source terms 


which are known as user-defined functions (UDF). The 
correctness of CFD model is very sensitive to the source 
term which can disturb the flow position of neighboring 
cells in both feed and permeate domains. 


Simulation setup and discretization 


Table 1 shows values of different properties used in 
this research work. Figure 3 illustrate the solution proce- 
dure for steady and unsteady (transient) problems. A 
UDF written in C programming language is developed 
and compiled with FVM solver. The library functions 
based on UDF build and load during the execution of a 
CFD model. Heat transfer from a hot air to a cellulose 
acetate based membrane is a problem that is aligned with 
the grid, the numerical diffusion will be certainly low, 
and so we can normally implement the first-order discre- 
tization scheme instead of the second-order scheme 
without any substantial loss of accuracy. The target of all 
discretization techniques like FVM is to develop a ma- 
thematical model to convert each term into an algebraic 
equation. Once implemented to complete control vo- 
lumes in a particular mesh, we attain a full linear system 
of equations that requires being solved. There are no 
common criteria for examining convergence. Definitions 
of residual that are advantageous for one type of problem 
are occasionally deceptive for other types of problems. 
So it is a good idea to examine convergence by judging 
residual stages besides observing related integrated quan- 
tities such as heat transfer coefficient. Intel Core™ i5 
with 4GB RAM is used to solve the steady and unsteady 
problem. For steady state problem, the CFD model is run 
for 100 iterations. For unsteady problem, the step size is 
0.1 seconds, a number of time steps are 200 and maxi- 
mum iterations per time step are 40. 


Results and discussion 


The default convergence criterion for most problems 
in FVM solver (ANSYS FLUENT™) is satisfactory. The 


Table 1 Values of different properties used in CFD modeling 


scaled residuals should decrease up to 10° for all equa- 
tions except the heat equation for which the scaled resi- 
duals should decrease up to 10°. Figure 4 shows that in 
this simulation the scaled residuals are decreased up to 
10°. Figure 5 illustrates the sum of facet values on a 
middle point. The middle point is created at the position 
(0.025, 0.05) in a rectangular grid of the membrane. 
ANSYS FLUENT™ is a cell centred solver, so it stores 
data based on the cell centered values. The only excep- 
tion is boundary conditions which are stored in the facet 
centers. As the value of temperature remains constant 
after 40 iterations (Figure 5), we can consider the solu- 
tion as converged. Figure 6 shows the contours of ther- 
mal diffusion from the feed side to a cellulose acetate 
based membrane due to a temperature gradient in a 
steady state problem. Figures 7-11 describe the contours 
of temperature diffusion in a membrane at various time 
steps for an unsteady state problem. The diffusion into 
the rectangular membrane from high to low temperature 
is clearly visible due to the use of mixed boundary condi- 
tion. This has been observed that maximum thermal dif- 
fusion occurs in first five seconds and the maximum sur- 
faced temperature of cellulose acetate based membrane 
reached up to 193 °C. The surface temperature increased 
up to 198 °C in twenty seconds. The process of heat 
transfer from air to the membrane is also due to the mo- 
tion of a fluid and diffusion of heat through the still 
boundary layer next to the membrane. 


Conclusion 


A CED solution for heat transfer problem was formu- 
lated to examine the performance of the thermal diffusion 
through the 2D rectangular cellulose acetate based mem- 
brane. The model was used to predict the steady and un- 
steady state problems in a membrane by solving the gen- 
eral scalar diffusion equation. In particular, we focused 
on high-temperature diffusion process in a membrane, 
which is industrially useful especially in chemical engi- 
neering. The FVM solver (ANSYS FLUENT™) pro- 


Properties Symbols Values Units References 
Thermal conductivity r 0.167 W/(m. K) [21] 
Specific heat of the material Cp 1451 J/(kg:K) [21] 
Heat transfer coefficient he W/(m?K) [22] 
Density (Air) p 1.225 kg/m? [23] 
Viscosity 4 1.7894x10° kg/m:s [23] 
°C 
Material temperature ġ K 
Surrounding fluid (Air) tem- °C 
j de [24] 
perature at feed side K 
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Fig. 3 Solution procedure for (a) steady state (b) unsteady 


state problems 
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Fig. 6 Contours of temperature for 
steady state problem in a mem- 


brane 
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Fig. 7 Contours of temperature at time= 
1 second 
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Fig.4 Residuals profile of UDS 
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Fig. 5 Convergence history of sum of facet values on a mid- 
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Fig. 8 Contours of temperature at time= 
5 second 
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Fig. 9 Contours of temperature at time= 
10 seconds 


vided the implementation of user-defined functions 
(UDFs) which have to be written in the C programming 
language. These UDFs were afterwards compiled for 
their use during solver execution for the solution of a 
problem. The residuals for the CFD model were de- 
creased beyond the standard level in order to achieve 
more reliable results. The results for the steady and un- 
steady problems showed that the developed CFD 
modeling approach is capable of predicting heat transfer 
in the polymeric membrane. This model can be extended 
to use for the modeling of complex gas mixtures as well 
as mass transfer problems to narrow the gap between 
laboratory and industrial scales. 


References 


[1] Versteeg, H.K., An introduction to computational fluid 
dynamics the finite volume method, 2/E. 1995: Pearson 
Education India. 

[2] Marriott, J., E. Sorensen, and I. Bogle, Detailed mathe- 
matical modelling of membrane modules. Computers & 
Chemical Engineering, 2001. 25(4): p. 693-700. 

[3] Wiley, D.E. and D.F. Fletcher, Computational fluid dy- 
namics modelling of flow and permeation for pres- 
sure-driven membrane processes. Desalination, 2002. 
145(1): p. 183-186. 

[4] Huang, L. and M.T. Morrissey, Finite element analysis as 
a tool for crossflow membrane filter simulation. Journal 
of membrane science, 1999. 155(1): p. 19-30. 

[5] Ghidossi, R., D. Veyret, and P. Moulin, Computational 
fluid dynamics applied to membranes: State of the art and 
opportunities. Chemical Engineering and Processing: 


1.98e+02 
1.90e+02 
1.82e+02 
1.73e+02 
1.65e+02 
1.56e+02 
1.48e+02 
1.39e+02 
1.31e+02 
1.23e+02 
1.14e+02 
1.06e+02 
9.73e+01 
8.89e+01 
8.05e+01 
7.21e+01 
6.37e+01 
5.53e+01 
4.68e+01 
3.84e+01 
3.00e+01 


Fig. 10 Contours of temperature at time= 
15 seconds 


[6] 


[11] 


[12] 


[13] 


[14] 


ChinaX iva (FRAT 


Fig. 11 Contours of temperature at time= 
20 seconds 


Process Intensification, 2006. 45(6): p. 437-454. 

Bao, L. and GG Lipscomb, Well-developed mass transfer 
in axial flows through randomly packed fiber bundles 
with constant wall flux. Chemical Engineering Science, 
2002. 57(1): p. 125-132. 

Lipscomb, GG. and S. Sonalkar, Sources of Nonideal 
Flow Distribution and Their Effect on the Performance of 
Hollow Fiber Gas Separation Modules. Separation & Pu- 
rification Reviews, 2005. 33(1): p. 41-76. 

Takaba, H. and S.-i. Nakao, Computational fluid dynam- 
ics study on concentration polarization in H< sub> 
2</sub>/CO separation membranes. Journal of membrane 
science, 2005. 249(1): p. 83-88. 

Schwinge, J., D. Wiley, and D. Fletcher, Simulation of 
the flow around spacer filaments between narrow channel 
walls. 1. Hydrodynamics. Industrial & engineering che- 
mistry research, 2002. 41(12): p. 2977-2987. 

Fletcher, D. and D. Wiley, A computational fluids dy- 
namics study of buoyancy effects in reverse osmosis. 
Journal of membrane science, 2004. 245(1): p. 175-181. 
Cao, Z., D. Wiley, and A. Fane, CFD simulations of 
net-type turbulence promoters in a narrow channel. Jour- 
nal of membrane science, 2001. 185(2): p. 157-176. 
Geraldes, V.t., V. Semiao, and M.N. de Pinho, Flow and 
mass transfer modelling of nanofiltration. Journal of 
membrane science, 2001. 191(1): p. 109-128. 

Koutsou, C., S. Yiantsios, and A. Karabelas, Numerical 
simulation of the flow in a plane channel containing a pe- 
riodic array of cylindrical turbulence promoters. Journal 
of membrane science, 2004. 231(1): p. 81-90. 

Koutsou, C., S. Yiantsios, and A. Karabelas, Direct nu- 


[15] 


[16] 


[17] 


[18] 


merical simulation of flow in spacer-filled channels: Ef- 
fect of spacer geometrical characteristics. Journal of 
membrane science, 2007. 291(1): p. 53-69. 

Koros, W. and G. Fleming, Membrane-based gas separa- 
tion. Journal of membrane science, 1993. 83(1): p. 1-80. 
McLellan, B., et al., Hydrogen production and utilisation 
opportunities for Australia. International journal of hy- 
drogen energy, 2005. 30(6): p. 669-679. 

Zhang, J., et al., Experimental and simulation studies on 
concentration polarization in H< sub> 2</sub> enrich- 
ment by highly permeable and selective Pd membranes. 
Journal of membrane science, 2006. 274(1): p. 83-91. 
Abdel-Jawad, M., et al., Flowfields on feed and permeate 
sides of tubular molecular sieving silica (MSS) mem- 
branes. Journal of membrane science, 2007. 299(1): p. 


229-235. 

Fluent, A., Fluent 6.3 User Guides: 8.6 User-Defined 
Scalar (UDS) Diffusivity. Fluent Inc, 2006: p. 50-58. 
Fluent, A., Fluent 6.3 UDF Manual. Fluent Inc, 2006. 
Acetate, P.C. Properties: Cellulose Acetate. 2014 
[cited Accessed June 2014]; Available from http://www. 
azom.com/properties.aspx ArticleID=1461. 

Ranges, H.t. Heat transfer coefficient 2014 [cited Accessed 
June 2014; Available from http://www.engineeringtoolbox. 
com/ convective-heat-transfer-d_430.html. 

Fluent, A., Fluent 6.3 Tutorial Guides. Fluent Inc, 2006. 
Arthanareeswaran, G., et al., Synthesis, characterization 
and thermal studies on cellulose acetate membranes with 
additive. European polymer journal, 2004. 40(9): p. 2153- 
2159. 


